Possible chiral phase transition in two-dimensional solid ^He 
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We study a spin system with two- and four-spin exchange interactions on the triangular lattice, 
which is a possible model for the nuclear magnetism of solid ^He layers. It is found that a novel spin 
structure with scalar chiral order appears if the four-spin interaction is dominant. Ground-state 
properties are studied using the spin-wave approximation. A phase transition concerning the scalar 
chirality occurs at a finite temperature, even though the dimensionality of the system is two and 
the interaction has isotropic spin symmetry. Critical properties of this transition are studied with 
Monte Carlo simulations in the classical limit. 
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Solid •^He films adsorbed on the surface of graphite 
are purely two-dimensional (2D) magnetic systems with 
isotropic interactions . Recently the solid phase of the 
adsorbed ^He layer at a low density has been attracting 
extensive interest, since it is a typical example of frus- 
trated quantum spin systems; nuclei of "^He form a spin- 
1/2 quantum antiferromagnet on the triangular lattice. 

In a double-layer solid '^He film, only the second layer 
is responsible for its magnetism at mK region. Specific- 
heat measurement of the second layer showed a peculiar 
behavior |^ and the effective exchange coupling obtained 
from the susceptibility was revealed to be antiferromag- 
netic [0 at a coverage where the layer just solidifies. It is 
believed that the second layer forms a triangular lattice. 
Elser proposed that the second layer may decompose into 
two sublattices, and atoms on only one sublattice which 
form a kagome net are responsible for the magnetism . 
Recent specific-heat data, however, showed that all spins 
on the triangular lattice contribute to the magnetism and 
they cannot be fully separated to the kagome lattice and 
the other QJ^. It is not yet clear whether the inequiva- 
lence of the two sublattices gives relevant effects on the 
translational symmetry of the magnetic interactions. An- 
other example of 2D solid '^He is realized in a monolayer 
of solid ^He adsorbed on graphite preplated with '*He ^] 
or HD §. 

It is well established that the multi-spin exchange in- 
teractions are important in the nuclear magnetism of 
solid ^He These multi-spin interactions originate 

from particle ring-exchange of "^He atoms. It was shown 
by Thouless that the whole Hamiltonian has the form 
H = X]n(~l)" S JnPm where Pn denotes cyclic permu- 
tation of n spins and the exchange constant J„ is always 
negative According to numerical calculations for the 
triangular lattice [|lO[|ll| ], two- and three-spin exchanges 
are large and, furthermore, four- and six-spin exchanges 
are not small. As pointed out in Ref. jl^, the multi- 
spin interactions can create frustration, and hence they 
should be properly considered in a theoretical treatment. 
In our previous study [|l3|, we found that the four-spin 



interaction can produce novel phases with four-sublattice 
structures. 

In this Letter, we reveal effects of the four-spin inter- 
action and predict a chiral phase transition in a certain 
parameter region where four-spin interactions are domi- 
nant. It is found that the ground state of a spin model 
for the solid ^He layer has a scalar chiral order. The spin- 
wave approximation shows that quantum fluctuations are 
not strong enough to destroy the chiral order at T = 0. 
We hence expect that the real "^He layer shows a chiral 
phase transition at a finite temperature. Critical prop- 
erties of this transition are also studied in the classical 
limit. The critical exponent a of the specific heat clearly 
deviates from that of the 2D Ising model. To our knowl- 
edge, the present model gives the first example of a chiral 
phase transition in a 2D realistic spin system with 5*0(3) 
symmetry. 

We consider a spin model with two-, three- and four- 
spin exchange interactions on the triangular lattice, and 
assume that the interactions have the same translational 
symmetry as that of the triangular lattice. Since the 
three-spin permutations can be transformed into two- 
spin exchanges, the Hamiltonian can be written as 



H 



(1) 



where J — J^^ J-if^, K = — J4/4(> 0), and cr, denote 
Pauli matrices. The first and the second summations 
run over all pairs of nearest neighbors and all minimum 
diamond clusters, respectively. The explicit form of hp 
for four sites (1,2,3,4) is 



hp= (Ti ■ (Tj + {(Tl ■ (T2){cr3 ■ (Ti) 

l<i<j<4 

+ {(Tl ■ (T4){(T2 ■ 0-3) - {(Tl ■ (T3){(T2 ' (Ti), (2) 



where (1,3) and (2,4) are diagonal bonds of the dia- 
mond. Bernu et al. first studied this model discussing 
the temperature dependence of the specific heat and the 
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FIG. 1. Spin vectors of four sublattices in the tetrahe- 
dral structure (left) and tlie configuration of the sublattices 
a, b, c, d on the triangular lattice (right). 

magnetic susceptibility |0 . Roger estimated density de- 
pendence of interactions using the WKB approximation 
[0. It shows that J3 is dominant in the high-density 
region and | J2 1 increases faster than | J3 1 as the density is 
lowered. Hence the value of J changes rapidly depending 
on the density and J may vanish at a low-density region. 

In the previous paper , we studied the ground state 
of the Hamiltonian (|^) for various J with the mean-field 
approximation. For J < — 8if , the ground state shows 
the perfect ferromagnetism. For —8K < J < —8K/3, 
ground states have non-trivial degeneracy. One of the 
ground states has a four-sublattice structure, in which 
three spins are up and the other is down. There ex- 
ist other kinds of ground states with longer periods. In 
—8K/3 < J < 25K/3, the ground state has a four- 
sublattice structure with zero magnetization, which we 
call as the tetrahedral structure (see Fig. 1). This spin 
configuration is indeed proved to be the exact ground 
state for the finite region —K/2 < J < 2K in the clas- 
sical limit. (See Ref. [^.) For 25K/3 < J, the ground 
state is the so-called 120° structure. Thus novel phases 
appear due to the four-spin interactions. Among these 
states, the tetrahedral structure has an interesting prop- 
erty; it has the scalar chiral order and hence the ground- 
state manifold has two- fold degeneracy. In this Letter, 
we employ the following order-parameter operator of the 
chirality @ 



E 



(Ti ■ {(Tj X (Tfc) 



(3) 



where the summation runs over all unit triangles and 
the indices {i,j,k) are chosen in clockwise turn for each 
triangle. 

To study the chiral order, we concentrate on the case 
only with the four-spin interaction, i.e., J = 0. This 
system has strong frustration due to the four-spin ex- 
change and hence quantum effects should be properly 
treated. We employ the spin-wave approximation to in- 
vestigate the strength of quantum fluctuations. Using 
the Holstein-Primakoff transformation, we expand the 
Hamiltonian in terms of Bose operators up to bilinear 



terms from the tetrahedral structure. After the Bogoli- 
ubov transformation of four kinds of bosons, we obtain 



H = - 



4 



(4) 



where the first summation runs over the reduced Bril- 
louin zone of the four-sublattice structure. The frequency 
cji(fc) is given by 



c.l(fc) = ^(4V-Sfc'-3C,2)l/^ 



(5) 
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The others are Lj2{k) = uji{k + Qi), cu^ik) — uJi{k + ^2) 
and uj4{k) = a;i(fc -I- -I- ^2)1 where and ^2 ^-re re- 
ciprocal lattice vectors of the four-sublattice structure. 
Three frequencies uJi{k) with i = 2,3,4 are gapless at 
fe = and the other one, ui^k), is massive. The ground- 
state energy per site is estimated as Eg = — 7.4615i^, 
whereas Eg = —17K/3 in the classical limit. The sub- 
lattice magnetization and the chirality are evaluated in 
the same way. The estimate of the sublattice magneti- 
zation per spin is rris/S = 0.5937, where S = 1/2. The 
deviation from the classical value, Atos/S" — 0.4063, is 
small compared with the value Anis/S = 0.522 of the 
pure Heisenberg antiferromagnet on the triangular lattice 
JlSf , which indicates that quantum fluctuations are not 
very strong and the ground state has the sublattice mag- 
netic order. The value of the chirality per unit triangle is 
k/2N = 0.2858, whereas k/2N = 4/3\/3 in the classical 
limit. The chirality has 37% of its classical value and 
this ratio is larger than the third power of the sublattice 
magnetization. This shows that the chiral order is more 
stable against quantum fluctuations than the sublattice 
magnetization. A similar tendency was observed in the 
vector chiral order of the Heisenberg and XY antiferro- 
magnets on the triangular lattice m&. The spin- wave 
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FIG. 2. Temperature dependence of the specific heat at 
J = 0. The figure around the peak is extended in the inset. 

analysis of the uuud and 120° structure ground states in- 
dicates that the stable region for the tetrahedral ground 
state will become narrower than the mean field result 

The greatest significance of the chiral order is that it 
can exist even at a finite temperature since it is stable 
against the spin-wave fluctuations with long wave-length. 
The spin-wave approximation indeed shows that the es- 
timate of the chiral order is nonvanishing at sufficiently 
low temperatures. We thus expect a phase transition 
to occur at a finite temperature, which is accompanied 
by ordering of the scalar chirality. We study the finite- 
temperature properties in the classical limit using Monte 
Carlo simulations. We believe that critical properties 
of the phase transition at finite temperatures are gov- 
erned by thermal fluctuations and hence the quantum 
effects do not change its universality class, though they 
can change values of Tc and the order parameter. We 
treat cr^ as a classical unit vector. Monte Carlo simula- 
tions are performed with the Metropolis algorithm. If a 
spin flip is rejected, we randomly rotate the spin about 
the local molecular field. We construct finite-size systems 
with a unit cluster which has 12 sites. The system size 
is TV = 12L2 with L = 4, 6, 8, 12, 16, 21 with a periodic- 
boundary condition. After discarding initial 6000-50000 
Monte Carlo steps per spin (MCS) for equilibration, sub- 
sequent 12 X 10^-5 X 10^ MCS are used to calculate the 
average. Further details will be reported in Ref. jl^ . 

The specific-heat data (Fig. 2) show a sharp peak 
around Tp — LOOK ~ 1.05K and this peak diverges 
as L increases. We verified that this phase transition is 
not of first order, by studying the cumulant of the fourth 
moment of the energy 1 — (i?^)/3(-E'^)^. The response 
function of the chirality, x = shows strong di- 

vergence at Tp and the long-range order of the chiral- 
ity; Vv*^/-^! appears below this temperature (Fig. 3). 



FIG. 3. Temperature dependence of the scalar chirality at 
J = 0. 

Of course, there is no sublattice magnetic order at finite 
temperatures even in the chiral ordered phase. Thus the 
divergence of the specific heat corresponds to the second- 
order phase transition from the disordered phase to the 
chiral ordered one. This is the first example of the scalar- 
chiral phase transition in the 2D realistic spin model with 
50(3) symmetry. We estimate the critical exponent a of 
the specific heat, C{T) - |r - Td"". The maximum 
values of the specific heat are plotted in Fig. 4 for var- 
ious sizes. The finite-size scaling reveals that the peak 
value behaves as C(Tp(L)) ~ 1°'/". The plot fits well by 
setting a/iy = 0.9(1). We hence obtain a = 0.62(5) and 
ly = 0.69(3) using the scaling relation a + 2iy ~ d. This 
is clearly different from the ^oij-divergence, i.e. a = 0, of 
the 2D Ising model. This result might be thought pecu- 
liar, since the chirality is an Ising-type variable and the 
chiral phase transition is expected to have the Ising uni- 
versality. There is however a similar example; the four- 
spin interaction added to the Ising model changes the 
exponent a continuously from | [T9[ | . The four-spin inter- 
action may have a universal effect to change the critical 
exponent. Other values are estimated as fi/v = 0.13(2), 
= 1.75(5) and Tc = 0.9935(10) from the finite-size 
scaling of the order parameter and the response function 
p^ . The critical indices may depend on the parameter 
J/K as they do in the 2D Ising model with four-spin 
interactions. Our estimates of the exponents (2 — a)/i^, 
l3/u, and "f/u appear to equal to 2, 1/8, and 7/4, re- 
spectively, and hence the present transition still belongs 
to the Ising universality-class in the sense of the weak 
universality [^ . 

Even for the vector chirality in the fully frustrated 2D 
XY model, there remain controversies whether the chiral 
phase transition has the Ising universality or not. Since 
Miyashita and Shiba estimated a — Q using Monte 
Carlo simulations, this problem has been studied repeat- 
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FIG. 4. Finite-size scaling plot of the peak values of the 
specific heat. 

edly |2^. Though the discussions have not yet settled 
down, all estimates of a are less than 0.4. Our estimate 
is clearly different from these values. 

We mention the possibility that the chiral phase tran- 
sition may be observed experimentally in the solid '^He 
films. The '^He films behave as ferromagnets in a high- 
density region [p] p^j2^ and, by lowering the coverage, 
the exchange coupling derived from the susceptibil- 
ity turns antiferromagnetic flpsf. On the other hand, 
the high-temperature series expansion shows that = 
— 2(J -I- &K) Hence the negativeness of indi- 

cates that J > and not necessary that J is anti- 

ferromagnetic. For all densities, no magnetic phase tran- 
sition has been observed at finite temperatures p^ , p5[ . 
Comparing these results with the phase diagram of the 
present spin model we find that the low-density 

■^He layer seems to correspond to the intermediate phase 
{—%K < J < —8K/3) between the ferromagnetic phase 
and the tetrahedral-structure phase. According to the 
WKB approximation the value of | J|(= | J2/2 - J3I) 
decays fast whereas K increases as the density is low- 
ered. We hence expect that, in a certain lower-density 
(but well solidified) region, the four-spin exchange inter- 
action is dominant and the chiral phase transition may 
occur. At this phase transition, one can observe sharp 
divergence of the specific heat, as we have shown in this 
Letter, and a cusp-like behavior in the magnetic suscep- 
tibility (see Ref. (isj). Recently, it was reported that, by 
preplating HD on graphite, the solid phase of a monolayer 
of "^He can be stabilized and then the ■^He film can first 
solidify around p = 0.0555A~^ [[7|p5|], which is the lowest 
density ever observed. (The second layer of double-layer 
■^He solidifies around p2 = 0.064A~^.) We expect that 
this material may be a plausible candidate for showing 
the chiral phase transition. 
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